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Using a sample of 771.6 x 10® T(4S') decays collected by the Belle experiment at the KEKB 
e^e“ collider, we observe for the first time the transition T(4S') —>■ r]hb{lP) with the branching 
fraction 'B[T(45') — >■ rjhb{lPy\ = (2.18 ± 0.11 ± 0.18) x 10“® and we measure the hb{lP) mass 
~ (9899.3 ± 0.4 ± 1.0) MeV/c^, corresponding to the hyperfine splitting AMhf(IJ’) = 
(0.6 ± 0.4 ± 1.0) MeV/c^. Using the transition hb{lP) —>■ -yribilS), we measure the rib{lS) mass 
.^rj 5 (is) = (9400.7 ± 1.7 ± 1.6) MeV/c^, corresponding to AMhf( 15) = (59.6 ± 1.7 ± 1.6) MeV/c^, 
the ribilS) width r^^(is) = (Sl® ± 5) MeV/c^ and the branching fraction B[/ii,(lP) —>■ 7 ? 7 i,(l<S')] = 
(56 ± 8 ± 4)%. 

PACS numbers: 14.40.Pq,112.38.Qk,12.38.Qk,12.39.Hg,13.20.Gd 


The bottomonium system, comprising bound states of 
b and b quarks, has been studied extensively in the past 
mm- The recent observations of unexpected hadronic 
transitions from the = 1 states above the BB 
meson threshold, T(45') and T(55'), to lower mass bot- 
tomonia have opened new pathways to the elusive spin- 
singlet states, the hb{nP) and r]b{nS) [3l |4j, and chal¬ 
lenged theoretical descriptions, showing a large violation 
of the selection rules that apply to transitions below the 
threshold. 


Hadronic transitions between the lowest mass quarko- 
nium levels can be described using the QCD multipole ex¬ 
pansion (ME) [5HTU] . In this approach, the heavy quarks 
emit two gluons that subsequently transform into light 
hadrons. The tttt and rj transitions between the vector 
states proceed via emission of ElEl and E1M2 gluons, re¬ 
spectively. Therefore, ij transitions are highly suppressed 
as they require a spin flip of the heavy quark mnn]. In- 
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deed, the ratio of branching fractions 
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vS 

tt-kS 


(n, m) 


'R[T{nS) ^ r]T{mS)] 
®[T(nS') ^7r+7r-T(mS')] 


is measured to be small for low-lying states: 1) = 

(1.64 ± 0.23) X 10-3 [TSHIS] and 3^^®s(3,l) < 2.3 x 
10-3 [n]. 

Above the BB threshold, BaBar observed the tran¬ 
sition T(45') —>■ 'qT{\S) with the unexpectedly large 
branching fraction of (1.96±0.28) x 10-“^, corresponding 
to 3illsi^,l) = 2.41 ±0.42 m- This apparent viola¬ 
tion of the heavy quark spin-symmetry was explained by 
the contribution of B meson loops or, equivalently, by 
the presence of a four-quark BB component inside the 
T(4S') wave function [T71[T8]. At the T(5S') energy, the 
anomaly is even more striking. The spin-flip processes 
T(5S') —>■ 7r7r/if,(lP, 2P) are found not to be suppressed 
with respect to the spin-symmetry preserving reactions 
T(5S') ^ 7r7rT(lS', 25') 0, and all the tttt transitions 
show the presence of new resonant structures [HI |20] 
that cannot be explained as conventional bottomonium 
states. 

Further insight into the mechanism of the hadronic 
transitions above the threshold can be gained by search¬ 
ing for the ElMl transition T(45) —)■ r]hi,{lP), which 
is predicted to have a branching fraction of the order of 

10-3 gl]. 

In this Letter, we report the first observation of the 
T(45) —>■ rjhi,{lP) transition and measurement of the 
ht,{lP) and rjb^lS) resonance parameters. Following the 
approach used for the observation of the hb{lP, 2P) pro¬ 
duction in e^'e- collisions at the T(55) energy 0 — by 
studying the inclusive tt+tt missing mass in hadronic 
events — we investigate the missing mass spectrum of 
r] mesons in the T(45) data sample. The missing mass 
is defined as Mmiss(f?) = \/{Pe+e- — P-qY^ where Pg+e- 
and Pri are the four-momenta of the colliding e“'"e- pair 
and the rj meson, respectively. 

The large sample of reconstructed /ib(lP) events allows 
us to measure its mass and, via the h},{lP) —>■ jrib^lS) 
transition, the mass and width of the r]b(lS). The latter 
are especially important since there is a 3.2 a discrepancy 
between the r]i,{lS) mass measurement by Belle using 
hb{lP,2P) —)■ 7?7b(15) transitions 0 and by BaBar and 
CLEO using T(25,35) —)■ 777t,(15) 

This analysis is based on the 711 fb-^ sample collected 
at the centre-of-mass energy of = 10.580 GeV/c^ by 
the Belle experiment [25l [26] at the KEKB asymmetric- 
energy e+e- collider p7H^ . corresponding to 771.6x 10® 
T(45) decays. Monte Carlo (MC) samples are generated 
using EvtGen m- The detector response is simulated 
with GEANT3 [31]. Separate MC samples are generated 
for each run period to account for the changing detector 
performance and accelerator conditions. 

Candidate events are requested to satisfy the stan¬ 
dard Belle hadronic selection [32], to have at least three 


charged tracks pointing towards the primary interaction 
vertex, a visible energy greater than 0.2-y/s, a total en¬ 
ergy deposition in the electromagnetic calorimeter (ECL) 
between O.l-^/s and 0.8y^, and a total momentum bal¬ 
anced along the z axis. Continuum e“'"e- —>■ qq events 
(where q G {u, d, s, c}) are suppressed by requiring P 2 , 
the ratio of the 2"'^ to 0*^^ Fox-Wolfram moment [35], to 
be less than 0.3. The 77 candidates are reconstructed in 
the dominant rj -G 77 channel. The 7 candidates are 
selected from energy deposits in the ECL that have a 
shape compatible with an electromagnetic shower, and 
are not associated with charged tracks. We investigate 
the absolute photon energy calibration using three cal¬ 
ibration samples: 7r° — >■ 77, rj -G 77, and D*^ -G P°7 
0. Comparing the peak position and the widths of the 
three calibration signals in the MC sample and in the 
data, as a function of the photon energy E, we deter¬ 
mine the photon energy correction Ten(^') and the reso¬ 
lution fudge factor 3fres{E). We observe 3fen{E) < 0.1% 
and 7res{E) ~ (±5 ± 3)% in the signal region, and ap¬ 
ply the corresponding correction to the MC samples. An 
energy threshold, ranging from 50 MeV to 95 MeV, is 
applied as a function of the polar angle to reject low en¬ 
ergy photons arising from the beam-related backgrounds. 
To reject photons from tt^ decays, 77 pairs having in¬ 
variant mass within 17 MeV/c^ of the nominal tt® mass 
[31] are identified as candidates and the corresponding 
photons are excluded from the rj reconstruction process. 
The angle 9 between the photon direction and that of 
the T(45) in the rj rest frame peaks at cos{9) « 1 for the 
remaining combinatorial background. We thus require 
cos(0) < 0.94 for the rj selection. All the selection crite¬ 
ria are optimized using the MC simulation by maximizing 
the figure of merit / = Ngig/^/Nsi^~+N^, where Agig 
and Abkg are the signal and background yields in the sig¬ 
nal region, respectively. The rj peak in the 77 invariant 
mass distribution, after the selection is applied, can be 
fit by a Crystal Ball (CB) [3S] probability density func¬ 
tion (PDF) with a resolution of 13 MeV/c^. Thus, 77 
pairs with an invariant mass within 26 MeV/of the 
nominal rj mass m,, |34| are selected as a signal sam¬ 
ple, while the candidates in the regions 39 MeV/c^ < 
\^{n) ~ ”^j)l < 52 MeV/c^ are used as control samples. 
To improve the M^issiv) resolution, a mass-constrained 
fit is performed on the rj candidates in both the signal 
and control regions. The resulting M^issiv) distribution 
is shown in the inset of Fig. The T(45) -G rjhi,{\P) 
and T(45) —>■ 77T(15) peaks in M^issiq) are modeled 
with a CB PDF, whose Gaussian core resolutions are 
fixed according to the MC simulation. The parameters 
of the non-Gaussian tails, which account for the effects 
of the soft Initial State Radiation (ISR), are calculated 
assuming the next-to-leading order formula for the ISR 
emission probability m and by modeling the T(45) as 
a Breit-Wigner resonance with F = (20.5 ± 2.5) MeV jr? 
[34] . The M^issirj) spectrum is fitted in two separate 
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FIG. 1. Mmiss(T/) distribution after the background subtraction. The solid blue curve shows the fit with the signal PDFs, 
while the dashed red curve represents the background only hypothesis. The inset shows the Mmiss(»?) distribution before the 
background subtraction. 


intervals: (9.30,9.70) GeV/c^ and (9.70,10.00) GeV/c^. 
In the first (second) interval, the combinatorial back¬ 
ground is described with a 6 **'-order (ll***) Ghebyshev 
polynomial. The polynomial order is determined maxi¬ 
mizing the credibility level of the fit and is validated using 
the sideband samples. Figure shows the background- 
subtracted M^issiv) distribution, with a bin size 50 times 
larger than that used for the fit. The credibility levels of 
the hts are 1% in the lower interval and 19% in the up¬ 
per one. The transition T(4S') —>• r]hb{lP) is observed 
with a statistical significance of Her, calculated using 
the prohle likelihood method [SS], and no signal is ob¬ 
served in the 77 -mass control regions. The hb{lP) yield 
is = 112469 ± 5537. From the position of the 

peak, we measure = (9899.3±0.4± 1.0) MeV/c^ 

(hereinafter the first error is statistical and the second is 
systematic). We calculate the branching fraction of the 
transition as 


®[T(4,5) ^ vhilP)] 


_ ^KjlP) _ 

^ri 4 S)<^rihbilP)'^[V 77 ] ’ 


where fV-f( 45 ) = (771.6 ± 10.6) x 10® is the number of 
T(4S'), SrjhbiiP) = (16.96 ± 1.12)% is the reconstruction 
efficiency and 23[77 — 77 ] = (39.41 ± 0.21)% [31]. We ob¬ 
tain ®[T(4S') ^ r]hb{lP)] = (2.18±0.11±0.18) x lO"®, in 
agreement with the available theoretical prediction [ 2 T] . 
No evidence of T(4S') — r]T{lS) is present, so we set 
the 90% Credibility Level (CL) upper limit 25[T(4S') —>■ 
r7T(15')] < 2.7 x 10“^, in agreement with the previous ex¬ 
perimental result by BaBar m- All the upper limits pre¬ 


sented in this work are obtained using the CLg technique 
[3a 00 ] and include systematic uncertainties. Using our 
measurement of we calculate the correspond¬ 

ing IP hyperfine splitting, defined as the difference be¬ 
tween the Xbj(lP) spin-averaged mass ^^^“^(ip) and the 
hb{lP) mass, and obtain AMhp(IP) = (-l-0.6±0.4±1.0) 
MeV/c^; the systematic error includes the uncertainty on 
the value of ^“^(ip) [ 2 l|- 

As validation of our measurement, we study the rj —>■ 
7 r“'" 7 r“ 7 r® mode. The tt® candidate is reconstructed from 
a 77 pair with invariant mass within 17 MeV/c^ of the 
nominal tt® mass [33] while the tt^ candidates tracks are 
required to be associated with the primary interaction 
vertex and not identified as kaons by the particle iden- 
tiheation algorithm. We observe an excess in the signal 
region with statistical signihcance of 3.5cr and measure 
®[T(4,5) ^ = (2.3 ± 0.6) x lO"®, 

which is in agreement with the result from the 77 mode. 

The contributions to the systematic uncertainty in our 
measurements are summarized in Table [l] To estimate 
them, we first vary — simultaneously — the fit ranges 
within ±100 MeV/c^ and the order of the background 
polynomial between 7 (4) and 14 ( 8 ) in the upper (lower) 
interval. The average variation of the htted parameters 
when the fitting conditions are so changed is adopted 
as the ht-range/model systematic uncertainty. Similarly, 
we vary the bin width between 0.1 and 1 MeV/c^ and 
we treat the corresponding average variations as the bin- 
width systematic error. The ISR modeling contribution 
is due to the T(4S') width uncertainty [33]. The pres- 
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TABLE I. Systematic uncertainties in the determination of 
‘B[T(4S') —>■ T]hh{lP)], in units of %, and on Mh^{ip), in units 
of MeV/c^ 


Source 



Fit range and backgronnd PDF order ±2.4 

± 0.1 

Bin width 

±2.5 

± 0.1 

ISR modeling 

± 2.8 

±0.7 

Peaking backgrounds 

±0.5 

±0.4 

7 energy calibration 

± 1.2 

±0.3 

Reconstruction efficiency 

± 6.6 

- 

Nr US) 

±1.4 

- 

Beam energy 

± 0.0 

±0.4 

B[?7 ^ 77 ] 

±0.5 

- 

Total 

± 8.2 

± 1.0 


ence of peaking backgrounds is studied using MC sam¬ 
ples of inclusive BB events and bottomonium transi¬ 
tions. While no peaking background due to B meson 
decay has been identified, the as-yet-unobserved transi¬ 
tions T(4S') — 77 T( 1 ^£)i^ 2 ) —>■ 77 ? 7 Tf(lS') can appear 
in the Mmiss(? 7 ) spectrum as a CB-shaped peaking struc¬ 
ture at Mmiss(? 7 ) = 9.877 GeV/c^ with resolution of 10.6 
MeV/c^. We take this effect into account by repeat¬ 
ing the fit with and without an additional CB compo¬ 
nent. No signal is observed and we obtain the upper 
limit on the product of branching fractions 1B[T(4S') —>■ 
77 T( 1374 i^ 2 )] X ®[T( 13 Di^ 2 ) ^ ?7T(15)] < 0.8 x lO'^ 
(90% CL). The uncertainty on the photon energy cal¬ 
ibration factors is determined by varying both J'eniE) 
and 3^res{E) within their errors. The uncertainty on the 
reconstruction efficiency includes contributions from sev¬ 
eral sources. Using 121.4 fb“^ collected at the T(5S') 
energy, the T(5S') —)■ 7 r+ 7 r“T( 2 S') transition is recon¬ 
structed; comparing the i ?2 shape obtained from this 
data sample with the simulation provides a ±3% uncer¬ 
tainty related to the continuum rejection. A ±1% uncer¬ 
tainty is assigned for the efficiency of the hadronic event 
selection. The uncertainty on the photon reconstruction 
efficiency is estimated using D events to be 

±2.8% per photon, corresponding to ±5.6% per rj. The 
number of T (45') mesons is measured with a relative un¬ 
certainty of ±1.4% from the number of hadronic events 
after the subtraction of the continuum contribution us¬ 
ing off-resonance data. The absolute value of accelera¬ 
tor beam energies are calibrated by fully reconstructed 
B mesons. We observe a ±0.4 MeV/c^ fluctuation of 
^ hb ( iP ) due to the uncertainty on the B meson mass [S3] 
and a negligible effect on the branching ratio measure¬ 
ment. Finally, we include an uncertainty in the branching 
fraction due to the uncertainty in 23 [77 —>• 77] [33] • 

The study of the 7?b(15) is performed by reconstruct¬ 
ing the transitions T(45) —>■ rjhjjilP) —>• ri^rn,{lS). To 
extract the signal, we measure the number of T(45) —>■ 
r]hf,{lP) events lV/j^(ip) as a function of the variable 
Aikfjniss — d7miss(^7) 21fniiss(^), where A7miss(^7) 
the missing mass of the 777 system. The signal transi¬ 


tion will produce a peak in at nirj^^is) —rnktiiP)- 

The radiative photon arising from the /ib(lP) decay is 
reconstructed with the same criteria used in the rj —> 77 
selection, and the hi,{lP) yield in each AM^iss bin is 
measured with the fitting procedure described above. 
To assure the convergence of the Mjnissiv) bt in each 
AMyniss interval, the hb{lP) mass is fixed to 9899.3 
MeV/c^, the range is reduced to (9.80,9.95) GeV/c^ 
and the order of the background PDF polynomial is 
decreased to seven. The /ib(lP) yield as function of 
AMmiss, shown in Fig. exhibits an excess at AM^niss = 
~ with a statistical significance of 9cr. 

The ?7&(15) peak is described by the convolution of a 



FIG. 2. AMmiss distribution. The blue solid line shows our 
best fit, while the red, dashed line represents the background 
component. 

double-sided GB PDF, whose parameters are fixed ac¬ 
cording to the MC simulation, and a non-relativistic 
Breit-Wigner PDF that accounts for the natural 77;,(15) 
width. The background is described by an exponential. 
We measure = (—498.6 ± 1.7 ± 1.2) 

MeV/c^, = ( 8^5 ± 5) MeV/c^ and the num¬ 
ber of T(45) —)■ rihb{lP) —)■ 77777 ^( 15 ) events = 

33116 ± 4741. The credibility level of the fit is 50%. We 
calculate the branching fraction of the radiative transi¬ 
tion as 

nhilP) ^ IVbilS)] = ^ 

where _ i ggy j- 0.053 is the ratio of the 

^^7^7 {,(15) 

reconstruction efficiencies for T(45) —>■ r]hb{lP) and 
T(45) —>■ rihbilP) —)■ 77777t,(15). We obtain 'B[hb{lP) 
777b(15)] = (56 ± 8 ± 4)%. To estimate the systematic 
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uncertainties reported in Table [Til adopt the protocols 
discussed earlier. Uncertainties related to the M„iiss(??) 


TABLE II. Systematic uncertainties in the determination of 
the rib(lS) mass and width, in units of MeV/c^ and on S = 
—>■ 'yr]b{lS) , in units of %. 


Source 

AMmiaa 

r»jbCis) 

2 

fit range 

±0.8 

±3.0 

± 2.8 

Afmiss(7) bin width 

±0.0 

±0.1 

± 0.0 

Mmiasiv) polynomial order 

±0.1 

±1.9 

± 1.6 


±0.0 

±0.8 

± 1.1 

AMmias fit range 

±0.0 

±0.7 

± 2.2 

AMniiaa bin width 

±0.8 

±2.8 

± 5.2 

7 energy calibration 

±0.5 

±0.3 

± 1.2 

Reconstruction efficiency ratio 

- 

- 

± 2.8 

Total 

±1.2 

±4.7 

± 7.2 

are determined by changing 

the fit 

range, 

the 


fit 

width, the background-polynomial order and the fixed 
values of used in the fits. Similarly, the un¬ 

certainties arising from the AMmiss fit are studied by 
repeating it with different ranges and binning. The cal¬ 
ibration uncertainty accounts for the errors on the pho¬ 
ton energy calibration factors. The uncertainty due to 
the ratio of the reconstruction efficiencies arises entirely 
from the single-photon reconstruction efficiency. The 
rjbi^S) annihilates into two gluons, while the hb{lP) an¬ 
nihilates predominantly into three gluons, but the MC 
simulation indicates no significant difference in the i ?2 
shape. Therefore, the continuum suppression cut does 
not contribute to the uncertainty arising from the recon¬ 
struction efficiency ratio. We calculate the rn,{lS) mass 
as + AMmiss = (9400.7 ± 1.7 ± 1.6) 

MeV/c^. Assuming tox(is) = (9460.30 ± 0.26) MeV/c^ 
[331 , we calculate AMhf(15') = (59.6±1.7±1.6) MeV/c^. 

A summary of the results presented in this work is 


shown in Table III We report the first observation of a 
single-meson transition from spin-triplet to spin-singlet 
bottomonium states, T(45') — rjhb{lP). This process 
is found to be the strongest known transition from the 
T(4S') meson to lower bottomonium states. A new mea¬ 
surement of the hb{lP) mass is presented. The cor¬ 
responding IP hyperfine splitting is compatible with 
zero, which can be interpreted as evidence of the ab¬ 
sence of sizable long range spin-spin interactions. Ex¬ 
ploiting the radiative transition hb{lP) —>■ 7776 ( 16 '), we 
present a new measurement of the mass difference be¬ 
tween the hi,{lP) and the rib{lS) and, assuming our 
measurement of we calculate Our re¬ 

sult is in agreement with the value obtained with the 
T(56) —?> 7r+7r“h6(lP) —>■ 7r+7r“7776(16) process [3| but 
exhibits a discrepancy with the Ml-based measurements 
[ 22 IT 24 ] . From the theoretical point of view, our result 
is in agreement with the predictions of many potential 
models and lattice calculations sn, including the re¬ 
cent lattice result in Ref. [32] . Our measurement of 
'B[hb{lP) —>■ 776(16)] agrees with the theoretical pre¬ 


dictions [331131] ■ All the direct measurements presented 
in this work are independent of the previous results re¬ 
ported by Belle [3|, which were obtained by reconstruct¬ 
ing different transitions and using a different data sample. 
Furthermore, all the results except for AMhf(16) and 
AMhf(IP) are obtained within the analysis described 
herein and are uncorrelated with the existing world av¬ 
erages. 


TABLE III. Summary of the results of the searches for 
T(45) nhb{lP) and hb{lP) 77ft(16). 

Observable Value 

‘B[T(46) ^ 7 /ib(lP)] (2.18 ±0.11 ±0.18) x 10“^ 

®[/i 6 (lP) ^ 776(15)] (56±8±4)% 

M 6 ^(ip) (9899.3 ± 0.4 ± 1.0) MeV/c^ 

- MhuiP) (-498.6 ± 1.7 ± 1.2) MeV/c" 
r,,(is) ( 8 l ®5 ± 5) MeV/c^ 

(9400.7 ± 1.7 ± 1.6) MeV/c^ 
AMhf(15) (±59.6 ± 1.7 ± 1.6) MeV/c^ 

AMhf(IP) (± 0.6 ± 0.4 ± 1 . 0 ) MeV/c^ 
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